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ABSTRACT 

Context. As is the case of several other Be stars, Achernar is surrounded by an envelope, recently detected by near-lR interferometry. 
Aims. We search for the signature of circumstellar emission at distances of a few stellar radii from Achernar, in the thermal IR domain. 
Meti:)Ods. We obtained interferometric observations on three VLTl baselines in the A' band (8-13 yum), using the MIDI instrument. 
Results. From the measured visibilities, we derive the angular extension and flux contribution of the A' band circumstellar emission 
in the polar direction of Achernar. The interferometrically resolved polar envelope contributes 13.4 ± 2.5 % of the photospheric flux 
in the A' band, with a full width at half maximum of 9.9 ± 2.3 mas (w 6^*). This flux contribution is in good agreement with the 
photometric IR excess of 10-20% measured by fitting the spectral energy distribution. Due to our limited azimuth coverage, we can 
only establish an upper limit of 5-10% for the equatorial envelope. We compare the observed properties of the envelope with an 
existing model of this star computed with the SIMECA code. 

Conclusions. The observed extended emission in the thermal IR along the polar direction of Achernar is well reproduced by the 
existing SIMECA model. Already detected at 2.2 /im, this polar envelope is most probably an observational signature of the fast wind 
ejected by the hot polar caps of the star. 

Key words. Stars: individual: Achernar; Stars: emission-line, Be; Methods: observational; Techniques: interferometric 



1. Introduction 

The southern Be star Achernar (aEri, HD 10144) has received 
much interest since its strongly distorted photosphere was re- 
solved by long-baseline interferometry (Domiciano de Souza et 
■ al. 120031 ), with major and minor axes of respectively = 2.13 + 
' 0.05 and 1.51 + 0.02milliarcseconds (Kervella & Domiciano de 
. Souza |2006l hereafter K06). Due to its extremely fast rotation 
(vsin/ « 250km. s ') and consequent flattening, the von Zeipel 
effect (von Zeipel 1924) causes the polar caps to be overheated: 
the polar effective temperature could be higher than 20 000 K, 
compared to < 10 000 K at the equator (Jackson et al. 120041 
Kanaan et al. 120081 hereafter Ka08). The high radiative pres- 
sure at the poles creates a fast polar wind that was detected in 
the near infrared by K06, where its flux reaches 4.7 + 0.3 % of 
the photosphere. In addition to this circumstellar envelope (here- 
after CSE), Kervella & Domiciano de Souza (2007) discovered 
a close-in companion of Achernar, of spectral type A1V-A3V 
(Kervella et al. I2008I I. In the present Letter, we report new in- 
terferometric observations of Achernar in the thermal infrared 
domain, using the VLTI/MIDI instrument. After a description of 
our measurements (Sect. |2]i, we derive the contribution and typ- 
ical angular scale of the polar CSE of Achernar using a simple 
Gaussian model and compare them to SIMECA model predic- 
tions (Sect.O. 
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2. Observations 

2.1. Instrumental setup and data processing 

MIDI (Leinert et al. l2003l Ratzka et aL l2007b is the mid-infrared 
two-telescope beam combiner of the Very Large Telescope 
Interferometer (VLTI; Glindemann et al. l2004l l. This instrument 
is a classical Michelson interferometer working in the astro- 
nomical band (7.6-13.3yum). For the reported observations 
of Achernar, we used a prism with a spectral resolution of 
R - A/AA ^ 30 to obtain spectrally dispersed fringes. During the 
observations, the secondary mirrors of the two Unit Telescopes 
were chopping with a frequency of 2 Hz to properly sample the 
sky background. Achernar was observed in 2006 and 2007, us- 
ing three 8.2 m telescope baselines (UT1-UT4, UT1-UT2, and 
UT3-UT4), and the SCIPHOT mode of MIDI. In this observing 
mode, the photometry of each telescope is recorded simultane- 
ously with the interferometric signals, allowing a more accurate 
calibration of the visibilities. The average dates of each obser- 
vation, with the corresponding projected baseline length B and 
position angle PA are given in Table [1] Each star or calibrator 
observation corresponds to a time on target recording interfer- 
ometric fringes of 3 min, followed by approximately 5 min of 
photometric calibrations. For the raw data processing, we used 
two different software packages: MIA developed at the Max- 
Planck-Institut fiir Astronomic and EWS developed at the Leiden 
Observatory (MIA+EW^ version 1.5.2) in order to extract the 
calibrated squared visibilities V^(A) (Chesneau 2007). We found 
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Table 1. Log of the observations of Achernar and its calibrator. 



# 


Baseline / Date 


UTC 


Target 


B (m) 


PA (°) 




UT1-UT4 










A 


2006-11-06 


01:59:46 


Q-Eri 


129.39 


48.89 


B 


2006-11-06 


02:26:23 


HD9362 


129.66 


54.55 


C 


2006-11-07 


00:35:03 


Q-Eri 


130.22 


32.23 


D 


2006-11-07 


01:05:39 


HD9362 


129.93 


40.06 




UT1-UT2 










E 


2006-11-06 


00:32:53 


Q-Eri 


52.38 


4.94 


F 


2006-11-06 


01:07:56 


HD9362 


54.64 


10.84 


G 


2006-11-06 


03:15:57 


ffEri 


49.58 


29.69 


H 


2006-11-06 


04:08:42 


HD9362 


50.53 


35.36 




UT3-UT4 










I 


2006-12-07 


02:48:10 


aEri 


62.46 


126.54 


J 


2006-12-07 


03:11:45 


HD9362 


61.76 


132.33 


K* 


2007-06-28 


07:15:29 


HD9362 


39.28 


50.04 


L* 


2007-06-28 


07:51:38 


cfEri 


46.10 


52.25 



* Due to a photometric instability, #K and #L were rejected. 



a good agreement between the results of the MIA and EWS 
packages within the error bars. In the following we will refer to 
the results obtained with the EWS package. Our calibrator, 6 Phe 
(HD 9362, G9III), was chosen in the Cohen et al. (1999 ) cata- 
logue of spectrophotometric standards for infrared wavelengths. 
(5Phe is located relatively close to Achernar on the sky (8.2°), 
and is of comparable brightness in the band (9.5 Jy at 12yum 
vs. «i 16 Jy for Achernar). It is almost unresolved by the inter- 
ferometer in the band, with 0ld - 2.24 + 0.02 mas (Borde 
et al. 2002). The calibrated squared visibilities of Achernar are 
listed in online Table|2] As MIDI operates in the diffraction lim- 
ited regime of the UTs, the effective field of view diameter is 
equal to ^ 0.26" at A - 10//m, much larger than the angular 
size of Achernar and its CSE. An imaging campaign in the 
band (Kervella & Domiciano de Souza 2007) uncovered a faint 
companion (Achernar B) located s; 0.3" away from the star 
Considering that the flux contribution of B is only * 2% of A 
in the band, we assume in the following that its impact on the 
MIDI visibilities is negligible compared to their accuracy, and 
refer to Achernar A simply as "Achernar". 

2.2. Visibilities 

The calibrated visibilities of Achernar are presented in Fig. [1] 
One problem we encountered is that the spectra for our two 
UT1-UT2 baseline observations #E and #G originally reached 
"unphysical" values larger than unity. Such a behavior was no- 
ticed by Chesneau (2007) , and is caused by an incorrect esti- 
mation of the K coefficients that characterize the internal pho- 
tometric transmission of the instrument. The resulting visibility 
spectrum is affected by a multiplicative bias larger than unity, 
uniform with wavelength. We could identify the cause for this 
behavior in the data sets #F and #H. The derived interferometric 
transfer function from these calibrator observations is sys- 
tematically lower than observation #B of the same calibrator, 
that was made (in time) in between observations #F and #H. The 
ratio r(#X) = r2(#X)/r2(#B) is respectively r(#F) = 0.86+0.07 
and 7(#H) = 0.82 + 0.03. As the seeing was stable over all our 
observations obtained on 2006- 11 -06 (0.6 - 0.9" in the visible), 
such a variation of cannot be explained by changing atmo- 
spheric conditions. The instrumental cause for this bias on the 
calibrator spectrum is not identified, but could be linked to 
an instrumental polarization problem. Our other observations are 



apparently free of such a "/f-induced" bias. The fact that we ob- 
tain nearly the same spectrum for two different nights on the 
UT1-UT4 baseline (Fig.[T] top) gives confidence in their calibra- 
tion. For this reason, we uniformly multiplied the spectra ob- 
tained on the UT1-UT2 baseline by the relevant y factors, whose 
uncertainties are reflected in the error bars (Fig.[T] middle). 

Achernar shows a squared visibility deficit on the longer 
UT1-UT4 baseline (Fig. [U top, B ^ 130 m), that is almost 
aligned with the rotation axis of the star Such a deficit (AV^ ^ 
20 - 30%) can be explained by the presence of a resolved CSE 
component, as discussed in Sect. [3] The visibilities on the UT3- 
UT4 baseline along the direction of the equator of the star show 
almost no resolution. However, the shorter projected baseline 
length of Si 62 m reduces the sensitivity of the interferometer 
to moderately extended emission. 

2.3. Spectropliotometry 

The absolutely calibrated spectrum of Achernar presented in 
Fig. |2] was obtained by dividing the average MIDI spectrum 
by the average spectrum of its calibrator HD 9362, and then 
multiplying the result by the template spectrum from Cohen 
et al. (I1999I ). The agreement with the IRAS spectrum (Volk & 
Cohen 119891 ) is satisfactory, although an excess can be noticed 
in the IRAS data between 8 and 10 /vm, compared to the MIDI 
and ISO spectra. This could be attributed to a different activity 
level of the star for these two observations. According to Fig. 14 
of Ka08, our MIDI observations were obtained in a state of in- 
creasing Ho- emission of Achernar, after a minimum occuring 
around 2000-2002. The IRAS observations took place in 1983, 
a year during which Achernar was in a decreasing, moderate ac- 
tivity state (Balona et al. [T987b . The ISO SWS spectrum (Sloan 
et al. 2003 ) was obtained in 1996, when the star was also in a 
decreasing activity state (Vinicius et al. 120061 ). 

The spectral energy distribution (SED) model presented in 
Fig.[3]was taken from the database of Castelli & Kurucz (120031) 
using an average T^ff of 15 000K and \ogg = 3.5 (Levenhagen 
& Leister 2006; Lovekin et al. 2006), for solar metallicity. The 
average angular diameter was set to ^ld - 1 -79 mas in order to 
match the observed broadband photometry in the V band, taken 
as a fiducial value. We chose this band as the contribution from 
the CSE is expected to be small in the visible. This value is very 
close to the arithmetic average of the polar and equatorial an- 
gular diameters (1.82 mas) measured by K06. We are aware that 
this SED model is not physically realistic (it ignores in particular 
the von Zeipel effect), but Lovekin et al. (2006) showed that the 
deviation from a rotating star SED is reasonably small. The pho- 
tometry was taken from Ducati (2002 ) for the U to N broadband 
photometry, Thompson et al. (1978) for the UV, ISO (Kessler 
et al. 12003] ). IRAS (IPAC [19861 ), and COBE/DIRBE (Smith et 
al. '2004') for the IR. We also used the VLT/VI SIR ph otometry 
obtained by Kervella & Domiciano de Souza ( 120071 ). as well 
as the average flux measured with MIDI between 8 and 12 /urn 
(Fig. |2]). As shown in Fig. [3] (bottom), an excess emission of 
ss 10-20% ofthephotospheric flux is present around /I = lOfim, 
but not at 25 and 60 jim (IRAS measurements). 

3. Extended emission 

To estimate the relative flux contribution and spatial extension 
of the CSE of Achernar, we use the simple model of a uniformly 
bright photosphere surrounded by a Gaussian CSE, with a full 
width at half maximum p and a flux contribution relatively to the 
photosphere a = /cse//*. In this model, we assume that the CSE 



p. Kervella et al.: Thermal infrared interferometry of Achemar with VLTI/MIDI 



3 



extension is independent of wavelength over the band. This 
type of model was already used by K06, and the interested reader 
is referred to these authors' Sect. 3.3 to 3.5 for details. In the 
present work, we restrict the fitting process to the polar direction 
of Achernar, i.e. to the UT1-UT4 and UT1-UT2 basehnes, that 
are approximately aligned (Fig.[T]|. We fix the photospheric polar 
angular size to dpo\ = 1.51 mas (K06). The adjusted model is: 



^model(P' ^pol' y) = [ j 



(1) 



where the photospheric (V*) and CSE (Vcse) visibilities are: 



V*(0pol,v) 



2yi(7r0poi v) 



Vcse(p, v) = exp 



(npyf 



4 ln2 



(2) 



where p is the FWHM of the CSE, and v - B/A the spatial 
frequency of the interferometric measurement. The result of this 
fit is shown in Fig. H] The derived polar CSE parameters are a 
FWHM of p = 9.9+2.3 mas and a flux contribution of o- = 13.4+ 
2.5 % relatively to the photosphere, on average over the band. 
The reduced of the fit is satisfactory at 1.6. The measured 
CSE flux is comparable to the excess measured photometrically 

of 10-20% (Sect. O- 

Along the equatorial direction, the visibilities obtained on 
the UT3-UT4 baseline (Fig.[T] bottom) are only marginally dif- 
ferent from unity. We can only set an upper limit of » 5% to the 
CSE flux contribution at 8 pm and w 10% at 12 yum in the equato- 
rial direction of the star, at spatial scales larger than 20 - 30 mas. 
The possibility still exists that a significant equatorial CSE con- 
tribution exists in this direction, but it should then be smaller 
than w 10/?*. A comparison of the photometric infrared excess 
measured in Sect. 12. 31 and the resolved polar CSE flux observed 
with MIDI also leaves little space for a bright equatorial com- 
ponent. We can roughly estimate the maximum equatorial CSE 
emission to ~ 5 - 10% at all spatial scales, in the band. 

For the present discussion, we choose to check the agreement 
of KaOS's Achernar model with the parameters derived from our 
Gaussian CSE model fitting, rather than directly fit the SIMECA 
model to the MIDI data. The reason for this indirect approach 
is the limited coverage in baseline orientation and length of the 
MIDI data that make the convergence of the fit difficult. The 
model of Achernar was computed by Ka08 using the SIMECA 
code (Stee & Bittar 2001) for their analysis of VLTI/VINCI data. 
SIMECA has been developed to model the environment of active 
hot stars. It computes line profiles, spectral energy distributions 
(SEDs), and intensity maps in lines and the continuum, which 
can be directly compared to spectroscopic, photometric and high 
angular resolution observations. The best model of KaOS for the 
epoch of VINCI observations (2002-2003) is that of a polar wind 
with an opening angle of about 20° in the K band (the model pa- 
rameters are listed in Table 1 of KaOS). From this same SIMECA 
model, we derive in the band a half-maximum CSE radius of 
6.8/?*, giving an angular FWHM of x 10 mas considering the 
X! 7/?0 polar radius of Achernar and its n = 22.68 + 0.57 paral- 
lax (ESA 1997 ). The associated A^ band flux contribution of the 
CSE is 1 1 % of the photosphere. The model extension and flux 
contribution are both in excellent agreement with the results of 
our MIDI observations (Sect.[3]l. 

The SIMECA model predicts CSE fluxes of w 2x and ^ 1 1 x 
the photosphere at 25 and 60 /im. Fig.[3]shows no such excess in 
the IRAS photometry. This discrepancy may be due to the low 
activity of Achernar during IRAS observations, or more proba- 
bly to inaccuracies of the model in this wavelength range. 



4. Conclusion 

From new interferometric observations in the thermal infrared 
domain (8-13 yum), we resolve an extended CSE along the polar 
direction of Achernar, whose total flux is 13.4 + 2.5 % of the 
photosphere, with a FWHM of 9.9 + 2.3 mas (« 6 /?*). This flux 
contribution is consistent with the photometric infrared excess of 
10-20% measured in the same wavelength domain, and with the 
predictions by the SIMECA model of Ka08. This convergence 
strengthens the plausibility of the presence of a fast polar wind 
ejected from the overheated polar caps of the star According 
to the scenario proposed by Ka08, Achernar is currently in a 
disk formation phase. Unfortunately, our limited coverage of the 
equatorial direction of Achernar with the MIDI data restricts our 
sensitivity. We could only set an upper flux limit of ^ 10% for an 
extended CSE in this direction. This non-detection is consistent 
with the K06 results in the K band. 
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Fig. 4. Observed values of Achernar and adjusted photosphere +Gaussian envelope model (dashed curve) along the polar direction 
(separate averages of UT1-UT4 and UT1-UT2 baselines), as a function of spatial frequency. The polar photospheric is shown as 
a dotted curve. 
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Fig. 1. Squared visibilities of Achernar. The error bars are the 
statistical dispersion of the measurements from the EWS soft- 
ware. The position angles of each projected baseline are given in 
degrees. The inset diagrams show each baseline in the {Au, Av) 
plane, and the polar direction of Achernar (dashed line). Thanks 
to the spectral coverage of MIDI, a relatively broad range of spa- 
tial frequencies is sampled simultaneously for each baseline (see 
Fig-©- The dotted line in the upperplot shows the polar func- 
tion of Achernar, for B = 130 m. The UT1-UT2 values were 
renormalized (see Sect. |2.2| for details). 
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Fig. 2. Absolutely calibrated MIDI spectrum of Achernar, using 
HD 9362 as a spectrophotometric standard star, with the spectra 
from IRAS LRS and ISO SWS (PWS processing) superposed. 
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Fig. 3. Spectral energy distribution of Achernar The open 
squares are measurements from the literature. The MIDI and 
VISIR points (Kervella & Domiciano de Souza l2007)) are shown 
as sohd dots. 
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Table 2. MIDI squared visibiUties of Achemar. The position angle is counted positively East of North {N = 0°,E = 90°). 



Baseline 
Proj.fi (m) 
Pos. Angle (°) 


UT1-UT4 
129.39 
48.89 


UT1-UT4 
130.22 
32.23 


UT1-UT2 
52.38 
4.94 


UT1-UT2 
49.58 
29.69 


UT3-UT4 
62.46 
126.54 


A(jj.m) 


VH^) ± a- 


y2(A) ± 0- 


VHX) ± cr 


V^iA) ± a 





7, 


.623 


0, 


.770 


+ 





.022 





.711 


+ 





.031 


0, 


.947 


+ 





.079 





.879 


+ 





.041 


0, 


.938 


+ 


0, 


.042 


7, 


.695 


0, 


.771 


± 


0, 


.022 


0, 


.693 


± 





.028 


0, 


.946 


± 


0, 


.079 


0, 


.890 


± 


0, 


.042 


0, 


.959 


± 


0, 


.042 


7, 


.767 


0, 


.760 


± 





.028 





.703 


± 





.031 


0, 


.948 


± 





.081 





.847 


± 





.037 


0, 


.991 


± 


0, 


.052 


7, 


.839 


0, 


.791 


+ 





.032 





.732 


+ 





.030 


0, 


.936 


+ 





.079 





.867 


+ 





.036 


0, 


.992 


+ 


0, 


.051 


7, 


.911 


0, 


.786 


+ 





.023 





.694 


+ 





.026 


0, 


.925 


+ 





.082 





.873 


+ 





.039 


1, 


.043 


+ 


0, 


.072 


7, 


.982 


0, 


.793 


+ 





.028 





.708 


+ 





.026 


0, 


.921 


+ 





.077 





.875 


+ 





.038 


1, 


.005 


+ 


0, 


.065 


8 


.053 


0, 


.812 


± 


0, 


.029 


0, 


.726 


± 





.025 


0, 


.908 


± 


0, 


.078 


0, 


.888 


± 


0, 


.038 


1, 


.011 


± 


0, 


.062 


8 


.123 


0, 


.802 


+ 





.028 





.725 


+ 





.024 


0, 


.906 


+ 





.076 





.894 


+ 





.039 


0, 


.999 


± 


0. 


.058 


8 


.194 


0, 


.810 


+ 





.024 





.729 


+ 





.025 


0, 


.901 


+ 





.076 





.875 


+ 





.037 


0, 


.988 


+ 


0, 


.046 


8, 


.263 


0, 


.809 


+ 





.025 





.737 


+ 





.025 


0, 


.901 


+ 





.077 





.881 


+ 





.037 


0, 


.971 


+ 


0, 


.045 


8, 


.333 


0, 


.812 


+ 





.024 





.754 


+ 





.025 


0, 


.901 


+ 





.077 





.870 


+ 





.036 


0, 


.969 


+ 


0, 


.045 


8 


.402 
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